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The local interplay between the ionic and electronic transport in NiO is explored using correlative imaging 
by first-order reversal curve measurements in current-voltage and electrochemical strain microscopy. 
Electronic current and electromechanical response are observed in reversible and electroforming regime. 
These studies provide insight into local mechanisms of electroresistive phenomena in NiO and establish 
universal method to study interplay between the ionic and electronic transport and electrochemical 
transformations in mixed electronic-ionic conductors. 

For over a century, electrochemical phenomena in solids were intensively explored in the context of energy 
storage and conversion devices 1 ' 2 , electrochromic systems 3 , corrosion, and metallurgy. Recently, much 
attention has been attracted to coupling between ionic and electronic transport phenomena, as exemplified 
by electroresistive materials and devices 4 . While the early applications of electroresistive systems date back for 
more than century (e.g. electronic coherer 5 ) and much work on fundamental mechanisms of electroresistance was 
performed in 1950ies, this field now sees a renaissance both due to emergence of new materials synthesis and 
nanofabrication technologies and recognized limitations of the complementary metal- oxide semiconductor 
(CMOS) technologies 6 ' 7 . The electroresistive systems are now actively explored in the context of non-volatile 
highly scalable information storage, memristive and neuromorphic electronics 7-9 . 

The processes involved in formation and operation of electroresistive oxide memories in oxides are very 
complex and involve multiple electronic and ionic transport steps and electrochemical transformations, often 
in nanometer volumes of active devices or conductive filaments and their junctions. As such, they sensitively 
depend on external stimuli, including bias levels, boundary conditions at the electrodes, initial microstructure, 
and forming currents, etc. 

The general sequence of events in electroforming of a material can be represented as following. For small 
external biases, only electronic conductivity is observed. For higher biases, the mobile ion subsystem is polarized, 
leading to accumulation/depletion of ions at electrodes. Associated changes in conductance lead to gradual 
development of in-plane instabilities that eventually gives rise to nucleation of new phases and formation of 
conductive filaments 71011 . In some materials, further bias increase (typically under galvanostatic control) leads to 
soft breakdown and formation of conductive filament (forming). In others, the reaction front is stable and bulk 
resistive switching is observed. In both cases, formation of gaseous products and oxygen adsorption on surfaces 
are observed, suggesting the role of gas-solid reactivity. Much effort has been directed towards controlling the 
electroresistive phenomena, e.g. through layers with given stoichiometry and device fabrication 7 ' 91012 " 14 . 

The design and optimization of electroresistive devices necessitates understanding of local ionic, electronic, 
and thermal phenomena over all stages of resistive switching, to determine the role of intrinsic defects and 
inhomogeneities on these processes. Recently, a number of studies of conductive filament formation in memris- 
tors was studied by (scanning) transmission electron and scanning probe microscopies (SPM) at nearly atomic 
resolution, providing the information on the structure and electronic properties of the formed filament 15 " 20 . In 
particular, atomically controlled electrochemical reactions including nucleation and subsequent growth have 
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been recently explored and provided information on the electrochem- 
ical kinetics of the conductive filament formation 18 " 21 . More recently, 
conductive filament formation has been studied on the nanoscale 
using current atomic force microscopy 22 . However, there is still lack 
of information on the electronic and ionic phenomena that precede 
the formation of conductive filaments, including ionic polarization, 
possible gas-solid reactions, and nucleation of conductive phase. 

Here we report the combined studies of electronic and ionic 
dynamics in NiO using correlative multimodal SPM. The electro- 
chemical strain microscopy (ESM) demonstrates the presence of 
strong hysteretic electromechanical response, indicative of bias- 
induced ionic transport. We further combine the first-order reversal 
curves (FORC) measurements in electrochemical strain micro- 
scopy 23 " 25 and in I-V measurements to decouple reversible and irre- 
versible processes. This correlative multimodal FORC IV and ESM 
imaging allows probing local mechanisms of electroresistive phe- 
nomena in NiO and establishing universal methodology to study 
the interplay between the ionic and electronic transport and electro- 
chemical transformations in mixed electronic-ionic conductors 26 . 

Results 

Electrochemical strain microscopy and spectroscopy of NiO. As a 

model system, we have chosen NiO. In the stoichiometric form, NiO 
is a classical cubic charge-transfer insulator. The changes in oxygen 
stoichiometry or doping by e.g. Li readily induce metal-insulator 
transition. The former can be controlled by electric bias, and 
consequently NiO is one of the most popular resistive switching 
materials 710,1114 ' 27 . Here, a 50 nm polycrystalline NiO!_ x thin film 
was prepared by sputtering method on SrRu0 3 /SrTi0 3 substrate 
to explore electronic and ionic dynamics using SPM. The details of 
the sample preparation and electrical properties can be found 
elsewhere 27 . 

The ionic dynamics in NiO is probed by electrochemical strain 
microscopy (ESM) 23 " 25,28,29 . In ESM, application of electric bias to a 
SPM tip leads to redistribution of mobile ions in the tip field 
and possible gas-solid reactions. The changes in molar volume and 



associated electrochemical strains lead to surface deformation that is 
detected by an SPM tip in the form of high-frequency (0.1-1 MHz) 
ESM signal. Indeed, strain effects induced by oxygen reactions were 
observed with atomic resolution using SPM 30 . In ESM voltage spec- 
troscopy, the response is measured as a function of slowly (~1 Hz) 
changing triangular voltage waveform, providing information on the 
bias-induced ionic processes below the tip. These imaging and spec- 
troscopic modes are similar to those adopted in piezoresponse force 
microscopy (PFM) 31,32 , with the difference being in the signal gen- 
eration mechanisms involved in imaging (piezoelectric response vs. 
bias-induced ionic motion) and spectroscopies (polarization switch- 
ing vs. electrochemical reactions). 

Shown in Figure 1(d) is surface topography of the NiO film, exhib- 
iting a clear grain structure. The ESM line images are shown in 
Figures l(a, b), illustrating the presence of ESM response of the 
pristine film. The locally varied ESM hysteresis loops can be observed 
as shown in Figure 1(c). The ESM hysteresis loops can be further 
analyzed over a dense grid of points 33 to extract the maps of effective 
switching parameters, as shown in Figures l(e, f). The spatial maps of 
switching parameters in Figures l(e, f) also show locally different 
electrochemical performance 24 . 

The observation of the strong electromechanical response in NiO 
and its hysteresis under external bias sweep strongly indicates the 
presence of ionic phenomena. Indeed, electromechanical response 
per se can originate from the piezoelectric effect, deformation poten- 
tial effects, ionic motion and associated Vegard strains, and flexo- 
electric effect 34 " 36 , as well as more subtle electromechanical couplings 
enabled by symmetry breaking at surfaces 37 . In cubic NiO, the bulk 
piezoelectric effect (and hence ferroelectricity) can be excluded. The 
deformation potential and flexoelectric effects describe fast res- 
ponses of the ionic sub -system to changes in the local electrochemical 
potential and field gradients respectively, and can contribute to the 
high-frequency ESM response. However, observation of ESM hys- 
teresis loop is indicative of either thermodynamic bistability (e.g. 
similar to that in ferroelectrics), or slow ionic or domain dynamics. 
In NiO, the latter can be attributed to the tip -induced ionic motion. 
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We note that while some of the ionically-induced electromehcnical 
responses can be non-trivial, e.g. due to the renormalization of the 
first term or Ginzburg- Landau free energy for polar mode or flex- 
oeffect, the controlling factor is still ionic dynamic 36 . 

Note that observation of cyclostationary ESM hysteresis loops 
suggest that the bias-induced ionic process is reversible or largely 
reversible (thus specifically excluding electroforming regions). In the 
NiO sample studied here, the ESM measurements are possible only 
in the +/ — 15 V bias interval. The application of higher biases leads 
to electroforming and dielectric breakdown, with associated large- 
scale changes in conductance and surface structure, as explored 
below. 

Probing ionic and electronic dynamics by first-order reversal 
curve spectroscopies. The ESM studies of NiO provide some 
insight into the local ionic phenomena in this material, as reflected 
in the spatial variability of the ESM response and hysteresis loops. 
However, comprehensive characterization of electroresistance me- 
chanism brings forward two challenges, namely (a) probing of both 
ionic and electronic degrees of freedom locally and (b) exploring the 
evolution of these responses during the irreversible processes. 
Current imaging or ESM imaging under a constant dc or ac 
condition does not provide any information on the evolution of 
these responses. Furthermore, a single hysteresis loop IV and ESM 
imaging can provide only limited information under a specific case. 
Thus, to address these challenges, we develop an approach based on 
correlative multimodal probing of ionic and electronic behaviors 
utilizing FORC method 38 " 41 applied to ESM and current-voltage (I- 
V) measurements. 

Macroscopic I-V curve measurement is a well-established appro- 
ach for characterization of electroresistive systems. For low voltages, 
application of bias to the tip does not alter the chemical state of the 
system. In this case, the I-V curve is non-hysteretic, i.e. forward 
and reverse curves coincide. For higher biases, the electronic flow 



through the system induces (reversible) changes in chemical or 
physical state of the material, including ionic polarization or local 
reduction. The relaxation processes are generally slow and are assoc- 
iated with significant kinetic or thermodynamic hysteresis. Corres- 
pondingly, the I-V curves become hysteretic in nature. Finally, 
formation of conductive filaments is associated with sharp changes 
in conductivity, constituting electroresistive switching. 

To explore these phenomena systematically, we implement the 
FORC I-V measurements in SPM spectroscopic imaging mode. 
Here, current collected from the bottom electrode is measured for 
the voltage profile shown in Figure 2(c). The area of the I-V hysteresis 
loop opening is probed as a function of the peak bias. For low biases, 
the I-V curves are non-hysteretic (or currents are below detection 
limit), whereas the onset of ionic transport is associated with a hys- 
teresis loop opening (see Fig. 2(a)). Correspondingly, the hysteresis 
loop opening vs. peak bias provides insight into the bias-induced 
electrochemical phenomena. 

Similar approach can be applied in ESM and PFM. For PFM, the 
FORC method allows local Preisach densities to be identified 42 ' 43 . For 
ESM, the FORC provides information on the development of ionic 
activity as a function of bias and allows separating of reaction and 
transport processes 44 . 

Several relevant comments regarding these methods should be 
made. At a single spatial location, FORC ESM and FORC I-V can 
both be applied to explore reversible and irreversible electrochemical 
processes. However, imaging on dense spatial grids {i.e. volume 
affected by bias is larger than pixel spacing) is possible only for 
reversible (or, more rigorously, cyclostationary) processes, while 
only sparse grids can be used for irreversible processes and generally 
spatially resolved imaging is unlikely 45 . That said, the onset of irre- 
versibility can be readily detected from images (streaks and abrupt 
changes of contrast, etc.). 

The high locality of electric field produced by SPM tip strongly 
affects the measurement mechanisms. Since bias is local, it suppresses 
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Figure 2 | Average FORC type (a) I-V and (b) ESM loops as a function of varying bias (c); (d) corresponding loop area averaged over 1600 spatial points. 
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(albeit not obviates) the filamentary instabilities usual in uniform 
electric fields (e.g. compare to ferroelectric polarization switching in 
PFM 46 " 51 ). Secondly, the ESM provides information on phenomena 
adjacent to the tip, i.e. within the signal generation volume that is 
commensurate with the tip-surface contact area 52 " 54 . While this pro- 
cesses can be affected by global conductivity of the film 55 , these effects 
are secondary and ESM signal can be observed even on purely ionic 
conductors 23 . Finally, the ESM signal is generally proportional to the 
ionic concentration. 

In comparison, I-V measurements measure the overall resistance 
between the tip and bottom electrode. While the spreading resistance 
of the tip -surface junction provides the dominant contribution to the 
signal for uniform materials, situation can be different for inhomo- 
geneous systems. Examples include the dislocation mediated con- 
ductance in STO where single defect establishes conductive path in 
the insulating matrix 56 . Opposite case is a grain boundary system, in 
which grain boundaries at large separation from the tip-surface junc- 
tion control resistance measured by SPM 57 " 62 . The current depends 
on the local ionic concentration; however, direct interpretation of 
spreading resistance in terms of carrier concentration is possible only 
for Ohmic conductance. For realistic materials, the I-V curves are 
highly non-linear 63 " 65 , precluding direct quantification in terms of 
relevant materials properties. 

Local ionic and electronic phenomena in NiO: reversible regime. 

The discussion above illustrates that the information provided by the 
FORC I-V and FORC ESM methods can be expected to be highly 
complementary. Here, we combine these two approaches in a 
correlative multimodal imaging mode, as shown in Figure 2. 
Figure 2(a) illustrates the evolution of the FORC I-V curves for 
small biases (in the reversible region). The corresponding ESM 
FORC data is shown in Figure 2(b). The peak bias dependence of 
the hysteresis loop area for I-V and ESM data is shown in Figure 2(d). 

The I-V curves are strongly non-linear, exhibiting shape char- 
acteristic for the mixed electronic-ionic conductors 63 " 65 . The hyster- 
esis loop opening in the I-V curve can then directly be attributed to 
the onset of ionic motion below the probe. The characteristic loop 



shape with regions of negative susceptibility at the turning points (on 
decreasing the bias hysteresis loop continues to increase) (see blue 
arrow in Figure 2(b)) strongly indicates the kinetic nature of the 
observed behavior. Remarkably, qualitatively similar behaviors are 
observed in the ESM data (Figure 2(d)). Here, for small bias windows 
the hysteresis loops are essentially closed and have characteristic 
lozenge shape. On increasing the bias, the loops evolve towards more 
rectangular shape; however the regions with negative susceptibility 
are also visible. Notably, the behavior of IV and ESM loop opening vs. 
bias is also very similar, hinting at the similarity of underpinning 
mechanism. This behavior is readily understandable in reversible 
regime, when vacancy accumulation below the tip provide linear 
correction both to electronic conduction and ESM signal. We also 
note that ESM loop opening occurs earlier than I-V loop opening. 
This originates from the fact that ionic motion can be detected by 
ESM even at lower bias which is not sufficient for current flowing 
through the films 66 . 

We further explore spatial localization of electrochemical and 
ionic activity on the NiO surface. In this, measurements such as those 
in Figure 2 are performed at each point of a dense spatial grid (40 X 
40 pixels with 25 nm pixel size). Note that the grid size is such that 
only several points fall within single grain. However, attempts at 
imaging with higher pixel densities typically lead to the irreversible 
surface degradation and were unsuccessful. 

Shown in Figures 3 (a, d) are the ESM hysteresis loop opening 
maps for the 8 V and 12 V bias windows (note that FORC data set 
by definition contains similar images for all bias windows with ideal 
spatial alignment). The presence of regions of higher- and lower ESM 
activity with the characteristic size of several tens of nm, i.e. com- 
parable to single grain size is clearly seen. The corresponding images 
for the FORC I-V measurements are shown in Figures 3(b, c, e, and 
f). Similarly to the ESM maps, clear localization of hysteresis res- 
ponses associated with the individual grains is observed. 

Of immediate interest is correlation between the ESM and I-V 
hysteretic behavior. In many cases, such correlation can be established 
based on the qualitative similarity of the images (large scale fea- 
tures, etc.). This is possible for the images in Figure 3; however, small 
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Figure 4 | Correlation between ESM loop area and I-V loop area spatial 
maps for (a) negative and (b) positive biases. 

feature size renders such comparison unreliable and it lacks quant- 
itative power. To explore this behavior systematically, we perform 
cross-correlation analysis of the images. Here, the cross-correlation 
coefficient for two images is defined as: 



Cor(A,B) = 



E w En(Amn-A)(B mn -B) 



(E m En(^--A) 2 )(E m En(^ 



-Bf) 



where images are represented as matrices A and B of equal size taken 
over identical spatial grid, A = mean (A) and B = mean (B). 

The cross-correlation between two images acquired at different 
biases establishes the measure of resemblance between the two, with 
a correlation coefficient of T meaning that features in both images 
are identical, '—1' being opposite (so that a minimum in one corre- 
sponds to a maximum in the other), and £ 0' meaning that the images 
are unrelated. Note that this analysis is limited in extracting the 
underlying relationships between data - for example the cross-cor- 
relation between topographic image and its derivative will be close to 
'0' (i.e. lack of correlation does not indicate lack of dependence); more 
complex multivariate statistical techniques can be used in these cases. 

Shown in Figure 4 is the bias -dependent correlation between the 
ESM and I-V hysteresis signals defined as Qj = Cor(ESM(Vj), 
IV ( Vj) ), where ESM( V,) and IV ( V,-) are the ESM and I-V hysteresis 



loop opening images at given biases. The correlation was calculated 
for negative and positive branches of the I-V loops separately. For 
negative biases, correlation matrix Qj is close zero for small biases, as 
expected for largely noise-dominated images. However, for large 
biases the correlation values are close to 1. This suggest that regions 
which has strongly hysteretic behavior in ESM also have strong 
hysteretic behavior in I-V curves, indicative of the common origin 
of the signals. For the positive biases, the behavior is remarkably 
different. Here, in the 6-8 V interval the Qj elements are negative, 
suggesting that the ESM activity is anti- correlated with the I-V. These 
features can be observed in the spatial maps of ESM loop area and I-V 
loop area in the orange boxes of Fig. 3. We note that some regions, i.e. 
grain boundaries, show opposite contrast dependent on the mag- 
nitude of biases. This means that these regions have dissimilar 
mechanisms of electrochemical response. 

Irreversible regime. On increasing the bias window, the material 
and tip become prone to irreversible changes of conductance 
(electroforming and tip degradation). The critical biases for these 
phenomena are position dependent, and in many cases their onset 
can be detected by spurious changes of ESM and I-V response during 
spectroscopic image acquisition. Note that once electroforming 
happened for the first time, the size of the affected area is often 
larger than the characteristic grid spacing, and hence all remaining 
measurements are taken on the electroformed region. Here we 
specifically aim to explore this regime to understand the changes 
in correlated IV and ESM behavior during electroforming. 

This behavior is illustrated in a set of FORC I-V and FORC ESM 
data in Figure 5 acquired sequentially for bias windows higher then 
stability region (previous section). For the first 10 points, the system 
remains in the pristine state, as shown in Figure 5. The hysteresis 
loop behavior is close to that in Figure 2, except that due to the higher 
bias the loops are more saturated. Again, the I-V and ESM loop 
opening exhibit similar dynamics. 

However, on continuous cycling the electroforming process 
occurs, with drastic changes in both I-V and ESM behavior. The 
I-V curve becomes strongly rectifying and the maximal current 
increases by ~2 orders of magnitude. At the same time, the ESM 
hysteresis loops exhibit strong vertical shift, adopt almost square 
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Figure 5 | Average FORC type (a, d) I-V and (b, e) ESM loops and (c, f) corresponding loop areas for (a-c) 10 points without any changes in surface state 
and (d-f) subsequent 1600 points including electroforming. 
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Figure 6 | Spatial maps of (a, d) negative I-V loop area, (b, e) positive I-V loop area, and (c, f) ESM loop area for (a-c) 2 V and (d-f) 14 V, respectively. 

Correlation between ESM loop area and I-V loop area for (g) negative and (h) positive biases. 



shape without the negative susceptibility regions, and effective res- 
ponse decreases by ~ 1 order of magnitude. 

The corresponding hysteresis loop area maps are shown in 
Figure 6. The clear spatially correlated features are visible in both 
ESM and I-V loop area maps. However, in this case the correlation 
between ESM and I-V behavior has almost disappeared. We attribute 
this behavior to the forming process in NiO. For intermediate biases, 
the oxygen vacancies accumulate below the tip, giving rise to large 
ESM responses and pronounced hysteresis in I-V and ESM loops. 
However, the system remains in the non-conductive state. On further 
bias increase, the nucleation of metallic phase occurs, with associated 
increase of electronic conductivity and decrease of ionic mobility. 

Discussion 

We have explored the interplay between electronic and ionic activity 
during the electroforming in NiO using newly introduced correlative 
multimodal imaging by FORC measurements in I-V and ESM. While 
conventional current or ESM imaging can provide only limited 
information on the electrochemical performance, the correlative 
multimodal imaging approach can provide electrochemical informa- 
tion of both ionic and electronic degrees of freedom locally. For small 
biases, both I-V and ESM provide the information on progressive 
vacancy accumulation under the probe with associated increase of 
zero-bias resistance and ESM loop opening as well as the contrast 
formation mechanisms. For large biases, the electroforming process 
occurs with a drastic increase of electronic conductivity and decrease 
of ESM responses, indicative of the conductive phase nucleation. 
These studies can be extended to other systems, most notably ferro- 
electric, tunneling barriers, fuel cells, etc., which allow local probing 
based on the electronic and electromechanical phenomena. 

Methods 

Materials. A 50 nm polycrystalline NiO!_ x thin film was prepared by sputtering 
method on SrRu0 3 /SrTi0 3 substrate (substrate temperature: 500°C, working 
pressure of Ar + 0 2 mixture gas: 1.5 mTorr, 0 2 ratio: 7%) to explore electronic and 
ionic dynamics using SPM. The surface roughness and grain size are about 8.852 and 
122.3 nm, respectively. The details of the sample preparation and electrical properties 
can be found elsewhere 27 . 

SPM measurements. Ambient SPM studies were performed with a commercial 
system (Asylum Research Cypher) additionally equipped with a current amplifier 
(FEMTO DLPCA-200) and LabView/Matlab based software. ESM and I-V 
measurements were simultaneously carried out with applying biases to Pt/Cr coated 



tips (the radius of the tip: <25 nm, Budget sensors Multi75E-G). The spatial 
resolution of the ESM and I-V maps is 25 nm. 
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